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ABSTRACT 

We present K-band polarimctric images of several massive young stellar objects at resolutions 
~ 0.1-0.5 arcsec. The polarization vectors around these sources are nearly centro-symmetric, 
indicating they are dominating the illumination of each field. Three out of the four sources 
show elongated low-polarization structures passing through the centers, suggesting the presence 
of polarization disks. These structures and their surrounding reflection nebulae make up bipolar 
outflow/disk systems, supporting the collapse/ accretion scenario as their low-mass siblings. In 
particular, S140 IRS1 show well defined outflow cavity walls and a polarization disk which matches 
the direction of previously observed equatorial disk wind, thus confirming the polarization disk 
is actually the circumstellar disk. To date, a dozen massive protostellar objects show evidence 
for the existence of disks; our work add additional samples around MYSOs equivalent to early 
B-type stars. 

Subject headings: stars:formation — stars:circumstellar matter — ISM:reflection nebulae 



Introduction 



With the wide availability of sub-arcsec observations, the presence of disks around massive young stellar 
objects (MYSOs) has become a hot topic recently. The reason for this, is the dispute about the way how 
massive stars are formed. Tha t massive stars are formed in a scaled-up version of how low-mass stars form 
(jShu. Adams, fc Lizano 1119871 ) would be a natural thought. However, it has been proposed that, if the mass 
of the central star exceeds 8 M Q , the tre mendous radiation pressure from th e central star will halt the mass 



accretion and keep it from growing (e.g.. lKahn 



1974 



Palla fc Stabler L99lh . Such a consideration leads to 



an entirely different mechanism, i. e., massive stars may be formed through mergers of lower-mass protostars 
(jBonnell. Bate fc Zinnecker II1998T ). On the other hand, a number of so lutions have been proposed to account 
for the radia tion pressur e problem, such as non-isotropic accretion ([Yorke fc S onnhaltcr 2002) , different 
dust opacity (jKahn 1119741) . and redirection of radiation by disk/outflow system ([Krumholz. McKee. fc Klein 
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2005). The presence of disk/outflow system around MYSOs will provide key evidence to evaluate these two 
kind of scenarios. 

Disks around MYSOs have been detected over a wide range of wavelength s with different tracers 
(jPatel et al. 1 120051 1 Jiang et al.l boosl ; ICesaroni et"al~l hood iBeltran et al~ll2006h . These results suggest 
that up to early B type stars can be formed through a similar route as solar type stars. In this work we 
report sub-arcsec near-infrared (NIR) polarimetric imaging observations which show three polarization disks 
that can be best interpreted as disks or toriods around the MYSOs. 



2. Observations and data reduction 



S140 IRS1 was observed on 2005 June 28, and the other thre e were observed on 2 006 November 5, 
all with the Coronagraphic Imager with Adaptive Optics (CIAO, Tamura et al. 2000f) mounted on the 



Subaru telescope. The polarization mode was set on (jTamura et al. 



2003h by placing a rotatable half-wave 



plate (images taken at waveplate angles of 0, 22.5, 45 and 67.5 degrees), upstream of the adaptive optics 
(AO) system, and a cooled wiregrid polarizer inside the CIAO cryostat. For each waveplate angle, the 
exposure times were 5 min for S140 IRS1, 4 min for S255 IRS1, 2 min for NGC7538 IRS1 and 6 min for 
IRAS23033+5951, respectively. The nights were clear and seeing was stable so that the point spread function 
of the optical system did not change much. For S140 IRS1, AO was not used due to a technical problem, 
resulting in a natural seeing, i.e., ~ 0.45" (FWHM); for the other targets, the AO system was applied, 
delivering a resolution of 0.15" to 0.2". The pixel size was set to 0.021" for all the observations. 



were calibrated with NGC 7538 IRS1 (|Yao fc Satolll9 99). Two schemes were applied to calculate the zero 



The images were reduced in the standard way with IRAF packa ges, as described in Jiang et al. I (|2005f ). 
The polarization data of S140 IRS1 were calibrated with A FGL 2591 (jTamura et al. Ill99ll) ; whilst the others 



angle correction of polarizatio n. One is to measure the p olarization of the po int sources of calibration, with 
the same aperture as done in iTamura et al. 1 ( 1991 ) and Yao fc Sato (1999). The other is to calculate the 
centroid of the polarization pattern of the assumed reflection areas. By adjusting the zero-angle correction, 
we can align the centroid to the center of the assumed illuminating source. In both schemes, the zero-angle 
correction agreed within one degree. Meanwhile, the assumption of centro-symmetric polarization pattern 
around NGC 7538 IRS1 can be tested by checking the non-calibrated polarization angles, which are roughly 
equal in every radiating line from the object. The polarization degrees were not calibrated because we 
used relative flux to calculate them. The instrumental polarization is negligible (less then 1%), and the 
polarization efficiency at K-band is about 97%. 



3. Results 

The main results are shown in Fig. 1. A common feature is noticed from inspection of the polarization 
vectors (left panels) for each field. In most of the fields, the polarization vectors are circling around the 
center, forming centro-symmetric patterns, suggesting the MYSOs are dominating the illumination of some 
or all of the fields. However, details of the polarization morphology are different from region to region, and 
are described below. 

S140 IRS1: The polarization vectors are perfectly centro-symmetric (Fig. la, left panel), suggesting 
that the source is dominating its surrounding area; no other sources affect the polarization. Two low- 
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polarization lanes pass through the central source. One is in the southwest-northeast at a position angle 
(p. a., measured east of north) ~ 45°, which is rather straight and narrow; the other is wider and a little 
curved to the north, and virtually perpendicular to the first one. In the pseudocolor image composed from 
intensity image and polarized flux image(Fig. la, right panel), two parabolically curved features, separated 
by the southwest-northeast lane, presumably representing reflection nebulosities, are seen opening towards 
the southeast and northwest, respectively. These two reflection nebulosities is probably the outflow cavity 
walls. A scattered photon is effectivel y polarized if the scatter ing angle (the angle between the incident and 
emerging direction) is near 90° (e.g., Code fc Whitney 199Sh . This is why we can only see the polarized 
light on two sides of the wall, while little polarization is observed in between. It is particularly true for large 
opening angles, as in the case of S140 IRS1. Th e wall profile to the northwest is more clearly observable , 



because it delineates a red-shifted outflow lobe ( Havashi et al. 1987 : Minchin. White &: Padman 1993), 



where extinction is larger so that we can only see the most efficiently polarized light on the two sides. The 
blue-shifted outflow to the southeast shows much brighter reflection nebulosity, consi stent with the result 



from the NIR speckle- interferogram images (jSchertl et al. 1 12000 : IWeigelt et al. 1 12002 ). The asymmetr y of 



the brightness between the nebulae, which had also been observed in CO lines (jHavashi et al. 1119870 . is 



probably a result of dust distribution around this object. The morphologies of the walls in dicate that the 



open ing angle is rather large (~150°), consistent with the poor outflow collimation ratio (jHavashi et al. 
Il987l ). 

The low-polarization lane in the southwest-northeast, therefore, re presents a polarizat ion disk (PD) 



Such a kind of feature is obs erved in intermediate-mass young sta rs ( Perrin et al. I I2004Q as well as in 



low- and high-mass protostars ([Lucas et al. I l2004r iJiang et al. Il2005h . iHoare I |2006h observed IRS1 in the 



centimeter continuum band, and found an elongated structure, which is interpreted as an equatorial disk 
wind. An over-plot of the ionized disk wind on the pseudocolor image shows that the directions of the 
PD and the disk wind matches each other (Fig. la, right panel). This evidence strongly suggests that the 
PD is actually a physical disk. The diameter of the PD is rather large (>2700 AU), implying the dusty 
disk is more extended than the ionized disk wind. Based on the bow-shock-like structures in the K'-band 



image. IWeigelt et al. I (|2002l ) proposed a quadrupolar outflow scenario for IRS1 (se e also lYao et al 



Wcigelt et al. 



1998|). 
2002), it 



However, since VLA4 is located along the path that the bow-shocks trace (Fig. 2 
is also possible that the 20° outflow is driven by VLA4. Given the configuration of the system we observed, 
a single bipolar outflow/disk system for IRS1 is more likely than multiple outflows. 

S255 IRS1: A bipolar nebula is located around S255 IRS1 (Fig. lb, right panel) at p. a. ^35°. High 
degrees of polarization (up to ^30%) in the nebula indicate that it should be a reflection nebula illuminated 
by the central source, IRS1, as illustrated by the polarization vectors (Fig. lb, left panel). I n the southwest 



part of the nebula there could be some contaminations-the outflow from IRS3 may overlap (jTamura et al 
19911 ). Still farther to the south (~ 3"from IRS1), a bright nebula is seen with scattering dominated by 



IRS3. In spite of this, the morphology of the bipolar nebula suggests an outflow fr om IRS1 with a sma ll 
opening angle (< 50°). Additional evidence is the H2 emission knots (S255:H2 2. iMiralles et al. 1119970 . 
which is located along the extension of the southern nebula (p. a. ~ 220 °), and which has been suggested 
to be excited by IRS1. The nebula is a little twisted in "S" shape, suggesting that it is a precessing outflow. 

A biconical low-polarization lane runs through IRS1, roughly perpendicular to the axis of the reflection 
nebula at p. a. ^110°, with a length of ^5700 AU. In Fig lb (right panel) we indicate the shape of the 
structure between two curves. The polarization vectors there are roughly parallel to the lane, but the whole 
patte rn of vectors is approxima tely elliptical. This structure is typic ally a PD, which has been well mod elled 
(e.g., Bastien &: Menard 1988h and observed in many cases (e.g., Whitney. Kenyon &: Gomez 1997 ). It 
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usually suggests a disk or a toriod around the central source. 



NGC7538 IRS1: Bright nebulosities are seen to the north and northeast of IRS1, but rather faint to 
the SE (Fig. lc, right panel). The polarization degrees to the north (up to 30%) are higher than to the 
south (~ 10%), and the polarized flux there show irregular structures. The polarization pattern of these 
nebulosities i s generally centro-sy mmetric with respect to IRS1 (Fig. lc, left panel). This confirms the 



suggestion bv lKraus et al. I (|2006f) that the bright patches to the northwest are illuminated by IRS1. Within 



~1" of IRS1, the polarization degrees are lower than the outside area, with a few low-polarization patches 
radiating out. There could be a PD surrounding this source, but we cannot tell which one it is from the 
complicated polarization morphology. 

From radio observations iGaume et al~ ( 1995h suggest a south- north outflow from IRS1. The reflection 
nebulosities to the nor thwest and north are therefore likely to be the cav ity wall of the blue-shifted outflow. 
Previous observations (jCampbell. fc Thompson 11984 1 iDavis et al. ~ll998h show the redshifted lobe is spatially 
coincident with the blue-shifted one which extends to the northwest, suggesting that the outflow is seen nearly 
polc-on. This would explain why we do not clearly see the southern counterpart of the reflection nebulosity 
where the extinction is much larger. 



From the methanol maser observation, iPestalozzi et al. I (|2004l ) suggested that there is a disk in the 
southeast-northwest direction. However, methanol masers can also be excited by high velocity gas. In light 
of our observation, it is also possible that the feature observed by Pestalozzi et al. represents o ne side of the 



cavit y walls. Such an interpretation does not conflict to the scenario of a precessmg outflow (jKraus et al. 
200fih . 



IRAS 23033+5951: This is the most complicated region of the four. In the pseudocolor image (Fig. 
Id, right panel), a conical nebulosity is seen, opening to the east. Several peaks are detected, but only the 
one marked by a plus (SI) is point-like (FWHM ~ 0.3"). To the east of SI, two other major peaks are 
detected. That they are not point-like and they have high polarizations suggests they are infrared reflection 
nebulosities (IRN). The northern one (IRN1) is mainly reflective; the polarization degrees are high (~20 %) 
and the vectors are facing towards SI, indicating it is illuminated by SI. The polarization degrees are lower 
(~10 %) in the southern one (IRN2), but the vectors also face to SI. The lower polarization could be caused 
either by a mixture of reflected light and self-radiation or by multiple illumination. To the northwest of SI, 
a faint reflection nebulosity is noticeable (IRN3), with relatively high polarization degrees (~30 %). The 
polarization vectors suggest this nebulosity is illuminated by another probably deeply embedded protostar 
that is not visible in the NIR. Since the property of this source is beyond the scope of this paper, a detailed 
analysis of it will be presented later (Fang et al., in preparation). 

Again, we see a low-polarization lane, which is more clearly seen in the polarization degree image (Fig Id, 
left panel), running roughly south-north (p. a. ^160°) and passing through the center of SI; the polarization 
vectors are generally aligned with the direction of the lane. As discussed above, this lane suggests the presence 
of a disk or a toroid. The length of this structure is about 6700 AU. The reflection nebulosities (IRN1 and 
IRN2) that are illuminated by SI, thus represent one direction of the outflow driv en by SI. This result is quite 



intrig uing since a bipolar outflow has been detected in the southeast-northwest (jKumar, Bachiller. fc Davis 



20021 ). Probably this outflow is driven by the deeply embedded protostar mentioned above. On the other 
hand, the outflow driven by SI has its own manifestation. A shocked H2 emission knot is found 4" west 
of SI (jKumar. Bachiller. fc Davis 1 12002). which is presumably excited by SI. It is interesting to note that 
the illuminating source, SI, is fainter than its illuminated reflection nebulosities in the K-band. A possible 
explanation is that the disk is seen edge-on so that the light directly from the central source is heavily 
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extincted while it can easily escape from the outflow cavity, then being reflected by the cavity walls. We 
note however, if only the extinction toward IRN1 and IRN2 is larger than that toward SI (A-vtiRNiy^-visi) 
>15), an observed effect can be reproduced. 



4. Implications and conclusions 

There are several possible causes that make a low-polarization lane passing through the bright source. 
In addition to a flattened structure around the illuminating source, a magnetic field may also produce a 
similar morphology by dichroic extinction. We checked the polarization of point sources in our observed 
fields, and find no large-scale alignment of polarization vectors in the direction of PDs. Another possibility 
is that, when a nearly spherical envelop illuminated by a source inside and a much brighter source outside, 
the envelop will show a bipolar polarization structure. However, this happens when the source outside is 
hundreds to thousands of times brighter than the source inside, so that the photons from the two sources 
shining on the envelop are roughly equal. It is not possible in our case since our sources are brighter than 10 3 
L/0 . It is also possible that a foreground filament which is just in front of the sources dilutes the polarization, 
but it is of very low probability and is in fact impossible to happen for all three targets. Therefore, taking 
the morphologies of the observed reflection nebulae into consideration, the most likely interpretation of the 
PDs here would be real disks or toriods. We note however, since the polarimetric observations do not give 
velocity information, it is not possible to distinguish between a Keplerian disk and a toriod structure. A 
confirmation of the conclusion have to await sensitive molecular line observations. 

Some parameters of the sources are listed in Table 1. Of the four sources, three (S140 IRS1, S255 IRS1 
and 123033) show evidence for the presence of a disk or toroid. All of them are deeply embedded (Ay > 15 for 
IRAS23033+5951, >20 for S255 IRS1 and >30 for S140 and NGC7538, Fig. 2a). They also show very large 
near infrared color excesses. Their near infrared colors suggest the evolutionary status could b e comparable 



to low-mass Class I objects. Except fo r S255 IRS1 , whos e spectral type is still controversial (jHever et al. 



1989c iHoward. Pipher. fcForrest 1 119971 : lltoh et al. I l200lh . these sources are MYSOs with masses ranging 



between 10 to 30 solar masses according to their spectral types. A J vs J-H plot would also suggest a similar 
conclusion about spectral types of these objects (Fig. 2b). 

Up to date, a dozen of MYSOs show evidence of the existence of disks ([Cesaroni et al. 200G; Z innecker fc Yorke 



2007). These candidate disks are found in various evolutionary status of their host objects, which are de- 



tected from (sub)millimeter wavelengths the NIR. However, the candidate disks h appen around MYSO s 



with their luminosities less that 10 5 L Q , typical of early B-type main-sequence stars ( Cesaroni et al. 2006), 
suggesting that disks might be common occurrence for this kind of MYSOs; we still lack of evidence whether 
still higher-mass YSOs host circumstcllar disks. Our observations add samples to this kind. And it is inter- 
esting to note that our most massive sample, NGC7538 IRS1 (~ 30 M Q ), does not show strong evidence of 
PD. 

We have also shown that high-resolution polarimetric imaging can provide a sensitive tool to detect 
circumstellar disks around MYSOs, even if the disks are too faint to be detectable directly. Our future work 
will concentrate on searching for polarization disks around MYSOs of higher masses and different evolutionary 
status. We aim to set up a sample for future observations with higher resolution and sensitivity. This kind 
of work will give crucial hints to answer the question at what mass and evolutionary phase the disks around 
MYSOs are truncated. 
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Table 1. Parameters of the sources 



Name 



D(kpc) 



L(L ) 



Sp type Sizc(AU) 



S140 IRS1 
S255 IRS1 
NGC 7538 IRS1 
IRAS 23033+5951 



O.f) 
2.4 
2.8 
3.5 



8.30E+04 
2.51E+04 



5.00E+03 



early B 
O,B0 
06 
B0. 5 



6700 



2700 
5700 



(1) 
(2,3) 
(4) 
(5) 



The projected length of the PD at adopted distance. 




2 Ref. Code 1. ILester et al. I |l98rt) : 2. iHoward. Pipher. feForrest I 



IWilliams. Fuller, fe Sridharan I 1120051 ) 
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\ 



,IRN3 



1RN2 



Fig. 1. — (left panels) Polarization degree images overlaid by polarization vectors (yellow dashes) and total 
intensity contour (blue curves). The angular and polarization scales are shown in the corners. The contours 
start from 10.93, 14.56, 11.34, 12.42 mag arcsec -2 , respectively, and decrease every 1.25 mag arcsec -2 . 
(right panels) Pseudocolor images composed of pure brightness images (red) and polarized brightness images 
(blue). A blue feature indicates the highly polarized nebulosity, which should be reflective. Red features are 
low polarization areas. The PD is indicated with double lines to delineate the edge of the disks (S140 IRS1 
and S255 IRS2), or a single line to indicate the position (IRAS23033+5951). For S140 IRS 1 The contour s 
of 5 GHz continuum emission representing the ionized equatorial disk wind are superimposed ( Hoare 20061 ). 
North is up and east is to the left. 



- 9 - 



3 



J/ /// S2*55 



X 

Hi 



NGC753 



140 



7 #3033 



1 2 
H-K 



4 

-2 B2( 



05 



N „5 •> 



68 



A2-. 



. 123033 



2 
4 

-05 0.0 



V -.NGC7538 
to ^ • 

\ v ;«.si4o 

v 



3.0 



Fig. 2. — a. Near infrared color-color diagram. The JHK magnitudes were obtained from the 2MASS 
database. Solid curves show the loci of main-sequence dwarfs and giants. The dotted lines indicate the 
reddening bands of main-sequence and T-Tauri stars, b. J vs J-H color-magnitude diagram of the sources. 
The solid curve indicates the zero-age main sequence stars with different spectral type. The dashed and 
dotted lines are drawn from B2 and 09 stars, parallel to the extinction vector, the scale of which is indicated 
by a cross on the dashed line. 



